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1 INTRODUCTION 



We present a new method to re-simulate the large scale 
structures of the Local Universe to a high precision. Re- 
simulating the dynamics of the large scale structures in 
our neighbourhood is an important step to understand 
the environmental effects that may affect galaxy forma- 
tion in our part of the Universe and to have better con- 
straints on the potential systematic effects in the obser - 
vati ons. Previous attem pts, such as iKolatt et ail (ll996l ) 
and iMathis et alj (|2002r ). used a global constraint of the 
smoothed density field, which is then transported back in 
time us ing Bernoulli-Zerdovich equation in a 80/t _1 Mpc 
volume. iNaravanan et all |200ll ) used PSCz and a pertur- 
bative expansion modelling of the evolution of the den- 
sity field to recover i nitial conditions in a 50/i _1 Mpc vol- 
ume. iFontanot et all (120031 ) used NOG and PSCz and a 
semi-Lagrangian reconstruction procedure to obtain the re- 
constructed displacements in a 6 0fc~ 1 Mpc vol u me. A dif- 
ferent approach was followed by iKlypin et all (|2003l ) and 
iKravtsov et~aH |2002i ) who used Mark III peculiar velocity 
survey and applied linear theory to obtain pure dynamical 
constraints, though rather sparse and in a 30/t _1 Mpc vol- 
ume. Related constrained simulations were then perform ed 
in the same manner (e.g iMartinez-Vaquero et al.l 120071 ) to 
study specific local environmental effects. All these recon- 
structions managed to capture the main large structures but 
the smaller scales Mpc) were not well recovered. 

In this work, we use the Monge- A mpere- K ant or ovit ch 
(MAK) orbit reconstruction method (Brenier et alj 120031 ; 
iMohavaee fc Tullvl 120051 ; lLavaux et ail |2008h to recover a 
dense set of constraints directly in Lagrangian coordinates. 
We use the catalogue o f the reconstructed orbits of the 
2MASS Redshift Su rvey (lHuchral|2000l; iHuchra et all 120051 ) 
already obtained in lLavaux et al. ( 2010l ). We put the con- 
straints object by object using a local adequate Lagrangian 
smoothing scale which depends on the estimated mass used 
for the reconstruction procedure. 



ABSTRACT 

We use the formalism of constrained Gaussian random fields to compute a precise 
large scale simulation of the 60/i _1 Mpc volume of our Local Universe. We derive the 
constraints from the reconstructed peculiar velocities of the 2MASS Redshift Survey. 
We obtain a correlation of 0.97 between the log-density field of the dark matter distri- 
bution of the simulation and the log-density of observed galaxies of the Local Universe. 
We achieve a good comparison of the simulated velocity field to the observed velocity 
field obtained from the galaxy distances of the NBG-3k. At the end, we compare the 
two-point correlation function of both the 2MRS galaxies and of the dark matter par- 
ticles of the simulation. We conclude that this method is a very promising technique 
for exploring the dynamics and structure of the Universe in our neighbourhood. 

In Section[2] we present the details of the used methods. 
In Section [31 we specify the set up of initial conditions from 
the 2MASS Galaxy Redshift survey. In Section [4] we discuss 
our results and in Section [5] we conclude. 



2 METHODS 

In this Section, we describe the general methods used to 
compute and apply the constraints on the initial conditions 
of the Af-body simulation. In Section [2.11 we recall the the- 
ory of constrained Gaussian realisation. In Section 12.21 we 
recall the fundamentals of the MAK reconstruction. 



2.1 Gaussian constrained realisation 

To build the initial conditions of our simulation, we 
are using a constrained Gaussian re alisation algorithm 
(|van de Weygaert fc Bertschingei] 1 19961). derived from the 
Hoffman-Ribak method (|Hoffman fc Ribaklll99ll ). This al- 
gorithm allows us to put any linear constraint of a Gaussian 
random density field, which corresponds to the initial con- 
dition of an A^-body simulation. We give here a summary 
of the steps of the algorithm. We first compute the variance 
matrix Vj j of the kernel of the constraints: 



r, 



d 3 kP(fc)A:*(k) J R" J (k) 



(1) 



with Kj(\s) the kernel window functi o n giv en in appendix 
F of Ivan de Weygaert fc Bertschingerl (|l996h . We build the 
Fourier transform <5 r (k) of a Gaussian random field and then 
compute the constrained realisation 



«5(k) = l r (k) + Y, Vrffa - J?) Ki (k) 



(2) 
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with Cj the value of the constraint, iV c the number of con- 
straints and 

cf = J d 3 kA'*(k)5 r (k). (3) 

After the inverse Fourier transform, the field <5(r) satisfies 
both the constraints and the power spectrum. Practically 
all the integrals are computed on the lattice of the finite 
discrete Fourier coefficients. 

The type and number of constraints is left free but must 
be linearly obtained from the density field S. For example, 
one may put constraints on the local shape of the density 
field and its peculiar velocity at the same time. In this case, 
we have thus 10 constraints for the local shape of the den- 
sity field, which corresponds to its value, its first and second 
derivatives. Additionally, we have 3 constraints for the pe- 
culiar velocity. Once the constraints have been put on the 
density, the rest of the density field is generated to match a 
given cosmological power spectrum such that the statistics of 
the final initial conditions are statistically correct. We note 
that this constrained realisation algorithm is fundamentally 
in Lagrangian coordinates to set up initial conditions for 
the dynamics of matter particles. This is particularly true 
for the choice of the size of the filter in the kernel Ki, which 
must be understood in terms of Lagrangian coordinates. 

We implemented such an algorithm in a MPI parallel 
environment Q This code is capable of handling an arbitrary 
number of Gaussian constraints of any size. The memory 
and processor load is shared evenly among all nodes of the 
cluster, making this code scalable on a MPI cluster of com- 
puter. 

2.2 The MAK reconstruction 

We propose to use the results obtained through the Monge- 
Ampere-Kantorovitch (MAK) reconstruction to find the 
constraints to put on the initial conditions of an Af-body 
simulation. The MAK reconstruction is a non-linear scheme 
of recovering the peculiar velocity field from an evolved 
non-linear density field. It assumes that the dynamics has 
not suffered shell cr ossing. In that case, it has been shown 
{Bremer et al.1 120031 ) that it is possible to uniquely recon- 
struct the trajectories of the particles from initial to present 
time. Practically, the reconstruction correspo nds to finding 
the m apping a which minimises the quantity {Bremer et al.l 
l2003h 

JV 

^ = E(^«- X 2 , (4) 

i = l 

where the {x^} are the Eulerian coordinates of equal mass 
particles composing the present day non-linear density field, 
{q^} are the homogeneously distributed Lagrangian coordi- 
nates of this same distribution. This method is thus purely 
Lagrangian. Additionally, full mass conservation is hard en- 
coded in the method, removing the singularities that are 
present in linear theory in the neighbourhood of large mass 
concentration. The assumption that no-shell crossing occurs 

1 This code is available publicly as ICgcn on 
http : //www. iap . f r/users/lavaux/icgen/icgen.php 



is of course limited to larger scales, typically a few mega- 
parsecs, where structures are still in the Lagrangian per- 
turbative regime, called the Zel'dovich approximation. This 
modelling of the dynamics is sufficient for our purpose which 
is to compute initial conditions compatible with observations 
on large scales. 

To handle redshift space distortions, we use the 
Zel'dovich approximation to predict the peculiar velocities 
of particles at the moment of the reconstruction. The mod- 
ified action has been shown to work well on large scales for 
the reconstruction of pec uliar velocities l)Mohavaee fc Tullvl 
120051 ; lLavaux et al.|[201oh . 

This method is particularly well suited to work in con- 
junction with the Hoffman-Ribak method as we derive all 
quantities in Lagrangian coordinates. 



3 COMPUTING INITIAL CONDITIONS 

In this section, we describe the methodology used to go from 
the original 2MRS data-set to the finished initial conditions 
of the A^-body simulation. Section 13.11 describes the data- 
set. Section [3.21 describes how the constraints were applied 
to the random Gaussian field. At the end, Section T3.3I gives 
the details of the setup of the initial conditions, in particular 
the choice of the background cosmology. 



3.1 Dataset 

To get better constraints for the simulation of the Local 
Universe, we have used the 2MASS redshift survey, which 
has the advantage of being full-sky, compared to other sur- 
veys, while being complete up to ~60-80/i~ 1 Mpc. To derive 
the mass distribution for the orbit reconstruction, we as- 
sume a constant Mass-to-Light ratio. This hypothesis has 
already given successful results on the comparison between 
reconst ructed and observed peculiar velocities using this cat- 
alogue {Lavaux et al.|[2010h . 

We use the mass s plitter algorithm given in Appendix A 
of|L avaux et alJ l)200ct ) . The algorithm looks for the optimal 
equal mass splitting given that we have a fixed number of 
mass elements and that the relative error between the split- 
ted and the true mass must be as low as possible. After the 
reconstruction of the displacements of the equal mass parti- 
cles, the global displacement of each of the groups of galaxies 
is computed by a simple average of the displacements of the 
particles representing this group. 

To mitigate boundary effects in the inner part of the re- 
construction volume, we use a 100/i _1 Mpc deep reconstruc- 
tion, padded with an homog e neous distribution of particles, 
as specified in lLavaux et al.l {200$), and have extracted the 
objects within 60/i _1 Mpc for the Local Group. 

We have not used all of the objects for constraining the 
Gaussian realisation as it would have yielded a 19913 x 19913 
matrix to compute and invert. This would be expensive to 
compute exactly, and additionally inversion would proba- 
bly need to be regularised. Instead, we limit ourselves to 
constrain the reconstructed peculiar velocities of the bigger 
objects. We now describe the methodology used to constrain 
the Gaussian realisation. 
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3.2 Constraining initial conditions 

We constrain the realisation within 60ft _1 Mpc of the centre 
of the box, which corresponds to our mock observer, using 
the reconstructed velocities obtained from the MAK recon- 
struction. This reconstruction has been shown to recover 
with precisio n the linear regime from an evolved non-l inear 
density field jMohavaee et al.ll2006l ; lLavaux et alJl200St ). We 
use the reconstructed 3D displacement field in Lagrangian 
coordinates to constrain the motion of matter. We choose 
an isotropic Gaussian smoothing window whose radius, Ri, 
is determined by the mass of the considered object i in the 
catalogue and is given by: 



r, = 4= ( m 



1/3 



(5) 



We add a l/y/E to match the Gaussian filter scale to a top 
hat filter scale to second order in Fourier space. The top hat 
filter is a good indicator of the filtering scale as the pecu- 
liar velocities are computed by doing a simple average of all 
displacements of the particles of a given group of galaxies 
when computing the total displacement of this group. We ne- 
glect the impact of the actual shape of the Lagrangian patch 
of the considered objects. To avoid having too many con- 
straints, we enforce that the constraint should at least have 
a Lagrangian radius Ri > 1.5ft Mpc. This corresponds to 
considering only the peculiar velocities of groups of a mass 
greater than ~ 4 x 10 1 Mq . This allows us to reduce the size 
of the matrix to 3942x3942. In doing so we make use of the 
sensitivity of the peculiar velocities to large scale power in 
order to reduce the computational complexity while keeping 
the results intact. Indeed, for a typical ACDM cosmology, 
the velocity field is correlated at ~ 90% for distances less 
than ~ 8ft ~ Mpc. Consequently, we do not add much in- 
formation by increasing the number density of constraints. 
Instead, we may even decrease the numerical stability of the 
method, which relies on a large matrix inversion. We discuss 
the impact of keeping only one fourth of these constraints at 
the end of Section[4]3] We expect that increasing the density 
of constraints should not change our results significantly. We 
include in Appendix[X]a comparison between the simulated 
velocity field and the original constraints. 



3.3 Initial condition setup 

We use a lBardeen et al.l l| 19861 ) power spectrum with n s = 1, 
n m = 0.258, Ov = 0.742, H = 72 km.s^.Mpc -1 , cr 8 = 0.77. 
The time complexity increases with the number of con- 
straints No as 0(Nc x N s ), with Ng the mesh size for gen- 



erating initial conditions. We choose a mesh of N g 



256' 1 



particles, sampling a box of side L = 200ft~ 1 Mpc for gen- 
erating the constrained initial conditions. This allows us to 
shape the white noise at a resolution of ~lft _1 Mpc while 
keeping the time complexity low. At the end, we constrain 
1,314 positions in a ~60ft~ Mpc volume in Lagrangian co- 
ordinates, which corresponds to 3,942 constraints. The com- 
putation of the correlation matrix of the constraints at this 
resolution took a little less than a week on 128 processors 
on the Mercury cluster at NCSA. 



4 RESULTS 

In this section, we detail the results of the constrained 
simulation. W e used the iV-body simulation code RAMSES 
l|Tevssierll2002h to evolve the initial condition that are spec- 
ified in Section f3. 31 The background metric is chosen to be 
a standard ACDM cosmology with parameters Sl m = 0.258, 
Q A = 0.742, H = 72 km.s _1 .Mpc _1 . We limit the number 
of refinement of the AMR grid to seven, this corresponds 
to soften the gravitational field on the scale of ~ 6ft _1 kpc, 
which is the size of one element of the mesh. We first look 
qualitatively at the sky-maps of the distribution of the dark 
matter. We then check the adequacy of the simulated ve- 
locity field to the observed velocity field within a volume of 
3,000 km.s - . Finally, we compare quantitatively the simu- 
lated density field to the observed galaxy density field: by di- 
rect comparison and through the two point correlation func- 
tion. 



4.1 Sky- maps 

In the left column of Figure [TJ we represent the position 
of galaxies by redshift slices on an Aitoff projection of the 
real sky. In the right column, we represent the projected 
mass density of the simulation in these same slices. The 
colour contrast is corrected such that the maximum value 
represented in these sky map is proportional to l/2(dj+i — 
di)(di + di+i) 2 . This ensures that clusters are not excessively 
dominating the map at larger distances compared to maps 
of the structures at smaller distances. 

The visual comparison of the structures seen in the 
2MASS Redshift Survey and in the simulation looks qualita- 
tively in agreement. We reproduce nearly all the filamentary 
structure of the Local Universe. We now look at the details 
of the structures in each of the slices. 

In the 10 — 30ft _1 Mpc slice, Virgo is at its right place 
near the North Galactic Pole. At the south pole, there is a 
significant dark matter concentration that is probably linked 
to a set of galaxies present in the 2MRS in the same slice. We 
correctly recover the Fornax cluster at (I, b) ~ (236, 53) with 
the Puppis cluster. We see that the Local Void area is totally 
em pty as it is expect ed from recent peculiar velocity analy- 
sis (|Tullv et alj|2008h . This void unfortunately corresponds 
to a region that is currently obscured by the galactic bulge 
of our galaxy. It is present in the reconstructed velocity field 
because we assume that the galaxies above and below the 
galactic bulge were a fair representation of what is happen- 
ing in the direction of the bulge. As the number of galaxies is 
very low, it happens that we have introduced a Local Void of 
the right size, as expected from peculiar velocities. This void 
now appears in the simulation. Additionally, we recover all 
the observed filaments. In the right panel, we see the first 
tip of the Pisces cluster at (I, b) ~ (150,-30). We note a 
filament that goes from this cluster to the Great Attractor 
region in the Hercules sector. This filament is lightly visible 
in the distribution of galaxies on the left panel. 

In the 30 — 40ft -1 Mpc slice, the galaxy distribution fea- 
tures are almost completely reproduced in the dark matter 
density field of the iV-body simulation. The first tip of the 
Pisces cluster appear in that slice, even though it should 
be mostly placed in the the 50ft -1 Mpc region. In the Great 
Attractor direction, (1, 6) ~ (300, 15), we see that the simu- 
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Figure 1. Observed vs simulated sky-maps - Left column: Aitoff projection of the galaxies of the Two-Micron-All-Sky Redshift survey 
in different redshift slices, as indicated above the panel. Right column: Aitoff projected mass density in different redshift slices of the 
simulation, centred on the mock observer at the centre of the simulation. The contrast of the different slices have been adjusted for the 
distance of each slices. 



lation follows very precisely the distribution of galaxies on 
the left. The two concentrations of galaxies in the middle of 
the sky projection correspond to the Hydra-Centaurus su- 
percluster, as highlighted in the right panel. The filaments 
are correctly visually interconnected according to the galaxy 
distribution. The Local Void is still present in this slice. 

So far, the number of galaxies introduced to fill the 
Zone of Avoidance was small. In the 40 — 50/i -1 Mpc slice, 
this number increases. We recover all marked structures of 
the right panel of this slice. We also recover the structures 
introduced by the filling of the Zone of Avoidance. For ex- 
ample, the "S" shape in the direction (l,b) ~ (300,0) is 
recovered in the distribution of dark matter. So, the con- 
strained initial condition are working even better than ex- 
pected in that respect, even though this shape may be an 
artificial construction due to the Zone-of- Avoidance filling 
using mirrored galaxies. 

The galaxy distribution in the 50 — 60h~ Mpc slice is 
totally reproduced in the dark matter density field of the 
simulation. 



4.2 Velocity fields 

In Figure [2] we compare quantitatively the simulated line- 
of-sight component velocity field to the observed line-of- 
sight component of the velocity field. For obtaining the ob- 
serv ed velocity field w e used the extended NBG-3k cata- 
log jTullv et al.l [20081 ). The distance were matched to the 



group of galaxies of the 2MASS Redsh ift Survey galaxies 
so that 726 groups receive a distance (I Crook et al.l 120071 ; 
lLavaux et al.l 20101 '). We did an equivalent treatment for 
the simulation. We first grouped the particles according 
to a Friend-Of-Friend algorithm with a linking length of 
0.2n~ 1//3 , with n the m ean particle density of the simulation 
l|Efstathiou et al1ll98ct ). We applied a threshold of having 
a minimum of 8 particles. This grouping helps at reducing 
Finger-Of-God phenomena while keeping other redshift dis- 
tortion effects. 

The two velocity fields are co mputed in redshift space 
to avoid systema tic effects (e.g. iLvnden-Bell et all 1 19881 ; 
iDekel et al1ll990h . Moreover, to avoid being contaminated 
by Finger-of-god effects, we first computed groups using a 
Friend-Of-Friend algorithm in the simulation (linking length 
equal to 0.2 the mean inter-particle length). We also used the 
distance to grouped galaxies obtained from the anisotropic 
Friend-Of-Friend that we used to make the groups of galaxies 
in the 2MASS Redshi ft Survey. This a n isotro p ic Friend-of- 
Friend is described in Huchra fc Gellerl (| 19821 ): I Crook et al.l 
\200H ); lLavaux et all fcoid ). 

The two velocity fields agree well with each other, al- 
though there is a large scatter. This scatter has at least 
three origins: the modelling error of the dynamics in the re- 
construction of the orbits of the galaxies of the 2MRS, the 
modelling error of the Zel'dovich transport when applying 
the constraints coming from reconstruction and the mea- 
surement errors on the observed velocity field. 
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Figure 1. continued 
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Figure 2. Simulated peculiar velocities along the line-of-sight vs. 
observed peculiar velocities - We represent here the direct com- 
parison of the simulated velocity field against the observed veloc- 
ity field (line-of-sight component) within the 30/i — x Mpc volume. 
We subtracted the 30/i -1 Mpc bulk flow in the two cases. The 
filled contours corresponds to 50% (black), 68% (dark grey), 95% 
(grey), 99% (light grey) of the 30/i~ x Mpc volume. Both veloc- 
ity fields were adaptively smoothed with a spline of radius in the 
range [2.5; f3]/i — 1 Mpc, with a median of 7h~ ^-Mpc. 
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Figure 3. Correlation of the density field of the simulation to the 
density field of the 2MRS catalogue - We have extracted a volume 
limited sample of 60/i -1 Mpc from the 2MRS and the equivalent 
volume from the simulation. The density fields are both smoothed 
to 3.5/i Mpc with a Gaussian kernel. 
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Figure 4. Two point correlation function of the galaxies and of 
the dark matter in the simulation - We have computed the two- 
point correlation function of the dark matter particles of the sim- 
ulation (red) and of the galaxies of the 2MASS Redshift Survey 
(black) for the volume limited sample of 60h~ iJVIpc of volume. 
Wc represented the ratio b 2 = §2MR,s/£simulation m the inset plot. 
The error bars are estimated using a Jackknife test on 1,000 sub- 
volumes with 1,000 different realisations of the sampling noise of 
the volume. 



4.3 Density field 

In Fig. [3] we show the correlation of the density field of the 
simulation and of the 2MRS galaxy number, smoothed to 
3.5/i _1 Mpc with a Gaussian kernel. We define the correla- 
tion coefficient as the ratio: 

<log(p sim /p s T m ) log(p ob s/Pobs)} 



V ((log(p s im/p S im)) 2 )((log(p obs /pobs)) 2 ) ' 



(6) 



with psim the smoothed density field of the simulation, p b s 
the smoothed galaxy number density field in a volume lim- 
ited sample 2MASS Redshift Survey. The correlation coeffi- 
cient of the logarithm of the density field is 0.97. This corre- 
lation is not perfect. On the high end, we see that the main 
alignment of the scatter distribution is not strictly following 
the red diagonal. This could be an effect of the galaxy bias 
at that scale. Reducing the number of constraints on the ini- 
tial condition from 3,942 to ~1,000, either by selecting more 
massive objects or by randomly picking up one constraint of 
four, leads to a decrease of the correlation of the two density 
fields to 0.92. Additionally, we see a significant increase in 
the scatter in the high density region, the low density being 
apparently unaffected. The scatter in the simulated velocity 
field remains unaffected by this change. 



and the actual two point function of the galaxies of our Uni- 
verse may allow to have a better understanding of the galaxy 
distribution bias by direct measurement. We propose here 
a simple test to estimate the scale dependent bias in our 
specific part of the Universe by taking the ratio of these two 

functions. 

We use the estimator o f lLandv fc Szalavl (|l993f ) to com- 
pute the correlation function within a 60/i _1 Mpc volume: 



£00 = 



DD(s) - 2DR(s) + RR(s) 
RR(s) 



(7) 



with s the amplitude of the redshift separation. We have 
built a volume limited sample of the 2MRS such that it is 
complete up to a distance of 60ft" 1 Mpc from the observer, 
which corresponds to keeping 5,472 galaxies. In this test, we 
are not sensitive to cosmic variance effects as we are look- 
ing at the same part of the Universe both in the simulation 
and in the observations. So, the error bars on our correlation 
function should only reflect the sampling errors. To estimate 
these errors, we have sub-divided the 60/i _1 Mpc volume into 
1,000 sub- volumes and performed a Jackknife test by com- 
puting the correlation function on the whole 60/i _1 Mpc ex- 
cept the concerned volume. For each of the computation of 
the correlation function, we have used a different realisation 
of the sampling noise used to compute RR and DR. This 
generated 1,000 correlation functions compatible with the 
data. 

We have computed the correlation function of the sim- 
ulation in exactly the same way. We have first extracted 
randomly 5,472 particles, and we have applied the same 
cuts in volume as for the catalogue, including the Zone-of- 
Avoidance effect. 

The results are given in Fig. 0] The central point of 
the correlation function corresponds to the average of all 
the obtained correlation function by the Jackknife test and 
the Poisson sampling test. The vertical error bars reflect 
the standard deviation of £(s) according to its average. The 
horizontal error bar gives the size of the bin. 

We see that the two correlation function are closely fol- 
lowing each other except on the smaller scales. This is ex- 
pected as we implicitly assumed in the reconstruction that 
galaxies are fair tracers on large scales. We note that, in 
our case, the value of the bias for the galaxies of the 2MRS 
seems to be b = 0.87 ± 0.03. However, this result is prob- 
ably systematically affected by the particular value of erg 
for the realisation of the large scale mode of the simula- 
tion, which are more compatible with as = 0.83 instead of 
as = 0.77. This can introduce a small bias of the order of 
0.93, which corresponds to what we observe here. Addition- 
ally, the fair comparison between the simulated velocity field 
and the observed velocity field (Fig. [2]) reinforces this result 
as the velocity field is sensitive to the total matter distri- 
bution, even though we used a constant mass-to-luminosity 
ratio to derive the mass of each galaxy to reconstruct their 
orbits. 



4.4 Two point correlation function 

Another quantitative comparison between observations and 
simulation is given by the two point correlation function. 
Comparing the two point functions given by the dark mat- 
ter distribution of the re-simulation of our Local Universe 



5 CONCLUSION 

We have implemented a new method to re-simulate the Lo- 
cal Universe with a standard Af-body simulation code from 
the 2MRS catalogue of galaxies. This method allows us to 
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reproduce in many details the Large-scale structures of the 
Local Universe up to 60/i _1 Mpc and to a precision of a few 
megaparsecs. This method is easily applicable to larger vol- 
umes and surveys. We illustrate the quality of the simulated 
volume by direct comparison within a 30/i _1 Mpc volume of 
the simulated velocity field to the observed velocity field. 
Furthermore, we highlight a direct application of this simu- 
lation by the joint measurement of the two point correlation 
function in the simulation and in the actual observations of 
the positions of the galaxies of the 2MRS. 

This method opens us the possibility of studying quan- 
titatively both larger and smaller scales for our Local Uni- 
verse. On small scales, such a simulation would allow to run 
a simulation of the Local Group including environmental ef- 
fects. On large scales, we would like to extend this work to 
the quantitative comparison between CMB data and Large 
Scale structure data, with the major advantage, compared 
to earlier attempts, of taking into account non-linear dy- 
namics and clustering. Along with semi-analytic modelling 
of galaxy formation, it is possible to build a non-statistical 
simulation of the galaxies in our sky, which would enable us 
to make a direct quantitative comparison with observations. 
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Figure Al. Simulated velocity field vs. original reconstructed 
velocity field - We represent here a quantitative comparison of the 
originally reconstructed velocity field, which was used to constrain 
the initial conditions, against the final simulated velocity field. 
We used the same colour convention and the same smoothing 
procedure as in Figure [2] 

APPENDIX A: COMPARISON OF THE 
RECONSTRUCTED TO THE SIMULATED 
VELOCITY FIELD 

We represent in Figure IaTI a comparison of the original re- 
constructed velocity field, used to derive the constraints on 
the initial conditions, to the simulated velocity field, ob- 
tained by evolving the initial condition using an TV-body 
simulation. We note that the two fields are fairly equal up 
to some scatter that is probably due to small scale cluster- 
ing, which is not modelled by the MAK reconstruction, and 
the imperfections of the Zel'dovich approximation, used to 
derive the velocity field constraints. 



